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Anderson localization is a universal quantum feature caused by destructive interference. On the other hand 
chiral symmetry is a key ingredient in different problems of theoretical physics: from nonperturbative QCD to 
highly doped semiconductors. We investigate the interplay of these two phenomena in the context of a three- 
dimensional disordered system. We show that chiral symmetry induces an Anderson transition (AT) in the 
region close to the band center. Typical properties at the AT such as multifractality and critical statistics are 
quantitatively affected by this additional symmetry. The origin of the AT has been traced back to the power-law 
decay of the eigenstates; this feature may also be relevant in systems without chiral symmetry. 
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The combination of global symmetries and dimensionality 
provides us with a useful classification scheme for the study of 
Anderson localization. 1 According to the one-parameter scal- 
ing theory (ST), 2 all eigenstates of a disordered system are ex- 
ponentially localized in fewer than two dimensions in the ther- 
modynamic limit and undergo a localization-delocalization 
transition [or Anderson transition (AT)] in higher dimensions, 
if the strength of disorder is sufficiently weak. However, de- 
viations from ST predictions have been found in systems with 
additional discrete symmetries. A typical example is a dis- 
ordered lattice without on-site disorder but random hopping 
amplitudes. 3 Formally, the system may be considered as two 
independent sublattices with nonzero matrix elements con- 
necting only sites of different sublattices. The corresponding 
Hamiltonian is invariant under a simultaneous sign flip of an 
arbitrary eigenvalue and the eigenstate components of one of 
the associated sublattices. As a consequence of this discrete 
symmetry, referred to as chiral symmetry, eigenvalues of the 
Hamiltonian come in pairs of ±q, that is, the spectrum is in- 
variant around the point e = 0. This symmetry is an important 
ingredient in different problems of theoretical physics: from 
the QCD Dirac operator 4 to bosons in random media 5 and the 
spectrum of a highly doped semiconductor. 6 

From the definition of chiral symmetry it is clear that one 
must distinguish between eigenvalues near zero (the origin) 
and distant from it (the bulk). In the bulk, spectral correlations 
in the metallic limit are described by the universal results of 
random matrix theory (RMT), usually referred to as Wigner- 
Dyson (WD) statistics. 7 Weak-localization corrections to 
these universal results are obtained by mapping the localiza- 
tion problem onto a supersymmetric nonlinear a model^ The 
transition to localization occurring for stronger disorder is be- 
yond the reach of current analytical techniques. Numerical 
results suggest that the AT is characterized by multifractal 9 
eigenfunctions, a scale-invariant spectrum, 10 and spectral cor- 
relations (usually referred to as "critical statistics" 11 ) differ- 
ent from Poisson and WD statistics. 1012 These properties are 
supposed to be universal, namely, they do not depend on the 
microscopic details of the Hamiltonian but only on the dimen- 



sionality and symmetries of the system. 

In contrast, much less is known about the region close to 
the origin where the chiral symmetry plays a crucial role. An 
exception is the metallic limit, where the effect of this symme- 
try on level statistics 13 has been investigated in detailed with 
the help of the powerful analytical techniques of RMT. Be- 
yond this region we are still far from even a qualitative un- 
derstanding of the interplay between disorder and chiral sym- 
metry. As a general rule, chiral symmetry tends to delocal- 
ize eigenstates close to the origin since weak-localization cor- 
rections vanish. 14 Indeed in the one-dimensional case it can 
be rigorously proved 3 that, in disagreement with the ST pre- 
diction, eigenstates at the band center are not localized ex- 
ponentially. In two dimensions, localization properties seem 
to be very sensitive to the microscopic form of the random 
potential. 15 For the specific case of a weak Gaussian disorder, 
it is well established that eigenvectors at the band center re- 
main delocalized. 16 In three-dimensions (3D), for weak disor- 
der and time-reversal invariance, the results of Refs fl7lfl8lfl9l 
suggest that chiral symmetry may induce power-law localiza- 
tion of the zero-energy eigenstates. 

Surprisingly, with the exception of certain one-dimensional 
systems with-long range disorder, 20 ' 21 the transition to local- 
ization in systems with chiral symmetry has hardly been stud- 
ied despite its potential applications to the description of the 
chiral phase transition in QCD, 22 or the metal-insulator tran- 
sition in highly doped semiconductors. The present Brief Re- 
port is a step in this direction. We investigate the interplay be- 
tween chiral symmetry and localization in a three-dimensional 
system with short-range disorder. Our main finding is that chi- 
ral symmetry induces an AT in the region close to the origin. 
Eigenstates are power-law localized with an exponent that in- 
creases as we move from the origin. The AT occurs when this 
exponent matches the dimensionality of the space. We recall 
that, since the AT is to a great extent universal, the validity 
of our results does not depend on the microscopic details of 
the model but only on the dimensionality of the space and the 
(chiral) symmetry of the Hamiltonian. 

We consider the following 3D tight-binding Hamiltonian 
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with only off-diagonal disorder: 

W = ^[tye ie «at 0j -+H.c] , (1) 

where the sum is restricted to nearest neighbors and the op- 
erator etj (a|) destroys (creates) an electron at the ith site of 
the 3D cubic lattice. We break time-reversal invariance by 
introducing a random magnetic flux described by the above 
Peierls phase e l8i: > with dij uniformly distributed in the inter- 
val [— 7r, 71"]. Although we mainly focus on the broken time- 
reversal case, we have also checked that our main conclusions 
are valid if time-reversal invariance is preserved. The hopping 
integrals tij are real random variables satisfying the probabil- 
ity distribution 

VQntij) = l/W for - W/2 < InUj < W/2 , (2) 

and zero otherwise. The above exponential distribution pro- 
vides an effective way of setting the strength of off-diagonal 
disorder.— 

In order to proceed, we compute eigenvalues and eigenvec- 
tors of the Hamiltonian (fl~|i for different volumes L 3 by using 
standard numerical diagonalization techniques. The number 
of disorder realizations for each volume L 3 ranges from 10 4 
(L = 10) to 500 (L = 20). Hard- wall boundary conditions are 
imposed in order to better examine the asymptotic decay of 
eigenstates in real space (see below). Since we are interested 
in the effect of chiral symmetry we mainly focus on a small 
energy window close to the origin e = 0. The eigenvalues thus 
obtained are appropriately unfolded, i.e., they were rescaled 
so that the spectral density on a spectral window comprising 
several level spacings is unity. 

Our first task is to look for a mobility edge in a small 
spectral window close to the origin. The chosen interval, 
e G (10~ 4 ,10~ 3 ) in our case, is not important provided that 
it is close to the origin. For other intervals one gets similar 
results but for a different critical disorder. We recall that in 
general the value W = W c for which the system undergoes 
an AT is quite sensitive to the details of the Hamiltonian and 
consequently is not universal. In order to locate the mobility 
edge we use the finite size scaling method. 10 First we evalu- 
ate a certain spectral correlator for different sizes and disorder 
strengths W. Then we locate the mobility edge by finding 
the disorder W c such that the spectral correlator analyzed be- 
comes size independent^ 

In our case we investigate the level spacing distribution 
P(s) [the probability of finding two neighboring eigenvalues 
at a distance Si — (e.; + i — e»)/A, with A being the local 
mean level spacing]. The scaling behavior of P(s) is exam- 
ined through the following function of its variance; 2 ^- 

r)(L, W) = [var(s) - var WD ]/[ var P - var W D] , (3) 

which describes the relative deviation of var(s) from the WD 
limit. In Eq. (f3]l var(s) = (s ) — (s) 2 , where (• • •} denotes 
spectral and ensemble averaging, and varwD = 0.286 and 
varp = 1 are the variances of WD and Poisson statistics, re- 
spectively. Hence r\ = 1(0) for an insulator (metal). Any 
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FIG. 1 : Scaling variable r\ as a function of disorder W for different 
volumes; only eigenvalues in the window 10 -4 < e < 10 -3 have 
been considered. The system undergoes an AT at W — W c ~ 18±1. 
In all cases the statistical error (not shown) is ~ 0.003. 



other intermediate value of r\ in the thermodynamic limit is an 
indication of a mobility edge. 

In Fig. 1 we plot the W dependence of 77 for different 
system sizes. It is clear that, in the window of energy stud- 
ied e G (10~ 4 , 10 -3 ), the mobility edge signaling an AT is 
located around W c ~ 18. We have also found that W c in- 
creases as the energy window of interest gets closer to the 
origin. For a weaker (stronger) disorder, 77 tends slowly to the 
WD (Poisson) results. This slow convergence to WD statis- 
tics (see Fig. 1) suggests that eigenstates may still have some 
kind of structure even on the delocalized side of the transition. 
The analysis of eigenfunctions will show that this is the case. 
We have checked that the use of scaling variables other than 
r]ML2£. do not alter the results. For the sake of completeness 
we repeat the calculation for the case of time-reversal invari- 
ance, 9ij = in Eq. (Q}. We have also observed a mobility 
edge with similar properties but at a slightly weaker disorder 
W c ~ 15, in agreement with previous results for the standard 
3D AT. 26 Finally, we remark that the study of r/(W, L) for dif- 
ferent disorders and system sizes is equivalent to the standard 
renormalization group analysis of the dimensionless conduc- 
tance g. For instance, according to the one-parameter scal- 
ing theory, g should be scale invariant at the AT. Any spectral 
correlator, rj in particular, is a function of g only and, con- 
sequently, it should be scale invariant at the AT as well. We 
have chosen r\ instead of g for numerical reasons. The former 
is easier to calculate and provides much more stable results. 

We now study whether level statistics around the mobility 
edge are compatible with those of an AT. Spectral correlations 
at the AT (critical statistics) combine typical properties of a 
metal (WD) with those of an insulator (Poisson). Thus level 
repulsion P(s) — > for s — > 0, typical of WD statistics is still 
present at the AT. However, a long-range correlator such as the 
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FIG. 2: The Z ^> 1 behavior of the number variance E 2 (Z) and P(s) 
[inset (b)] in the spectral window e £ (10 -4 , 10 -3 ) at the critical 
point W c ~ 18 for different volumes. £ 2 (7) is linear with a slope 
(solid line) x ~ 0.49 ± 0.01, P(s) ~ s for s < 1, and both are 
size independent. These features are typical signatures of critical 
statistics. For W = 0.2 <C W c [inset (a)] the number variance 
(dotted line) is well described by the prediction of chiral RMT (solid 
line). 



the number variance £ 2 (£) = ((N e - (A^)) 2 ) ~ x^or £ » 1 
(7V^ is the number of eigenvalues in an interval of length £) is 
asymptotically linear, 12 as for an insulator E 2 (£) ~ £, but with 
a slope % < 1 (~ 0.27 for 3D AT). 

As shown in Fig. 2 all these features are also present in our 
model. However, we have found that the slope of the number 
variance x = 0.49 ± 0.01 is roughly twice that of the standard 
3D AT. We refer to the analysis of eigenvectors (see below) 
for a qualitative explanation of this feature. 

We have verified that the number variance in the limit of 
weak disorder agrees with the prediction 13 of chiral RMT [see 
inset (a) Fig. 2]. For weak enough disorder results depend 
on whether the time-reversal invariance is broken [inset (a) in 
Fig. 2] or not (not shown). As a general rule, as disorder 
is increased the magnetic flux gets weaker, for W > 10 we 
could not observe any difference. Additionally, we have found 
that, if chiral symmetry is broken, by for instance adding a 
weak diagonal disorder, level statistics are very close to the 
Poisson limit typical of an insulator. 

To conclude, the analysis of level statistics shows that, as 
a consequence of the chiral symmetry of the model, there is 
an AT close to the band center with properties similar but not 
equal to those of the standard 3D AT. 

One of the signatures of an AT transition is the multifrac- 
tality of the eigenstates. An eigenfunction is said to be multi- 
fractal (for an alternative definition see Ref. 27) if the eigen- 
function moments P q = J d d r\ip(r)\ 2q oc L - - ^ 9 " 1 ) present 
anomalous scaling with respect to the sample size L, where 
D q is a set of different exponents describing the transitioni 2 ^ 2 ^ 

After performing ensemble and spectral averaging in the 



spectral window e £ (10 4 , 10 3 ) we have found that the 
eigenfunctions are indeed multifractal (D^ = 0.69 ± 0.02, 
D 3 = 0.45 ± 0.01, and D 4 = 0.340 ±0.007). These values of 
D q are, roughly speaking, a factor of 2 smaller than the ones 
at the standard 3D AT. This is consistent with our results from 
level statistics, where a similar factor was found for the slope 
X of the number variance S 2 (£) ~ ■ Qualitatively, this nu- 
merical difference can be traced back to the chiral symmetry 
of our model: according to the ST, parameters such as x or 
D q are functions of the dimensionless conductance g c at the 
AT. Although an explicit expression is not known, it is ex- 
pected that, at least for sufficiently high dimensions, both D-2 
and 1 — x w iH t> e proportional to g c . On the other hand, it 
is well established 1416 that weak-localization corrections van- 
ish in systems with chiral symmetry. Hence, a stronger dis- 
order is needed to reach the AT region. As a consequence, 
the chiral g c will be smaller, in qualitative agreement with our 
results. Having discussed the details of the AT close to the 
origin, we now clarify the physical mechanism causing this 
transition. Chiral systems do not have weak-localization cor- 
rections, so the transition to localization must be in part due to 
some nonperturbative effect, such as tunneling combined with 
the progressive weakening of the effect of the chiral symmetry 
as we move away from the origin. On the other hand, the slow 
rate of convergence toward WD statistics observed in Fig. 1 
suggests that, although eigenstates seem to be delocalized in 
the thermodynamic limit, they still have some kind of local- 
ization center. This is in agreement with the results of Ref. 
IT7I where it was found that the zero-energy wave function 
of a time-reversed, weakly disordered 3D chiral system was 
power-law localized ijjjr) ~ \/r a ^ w \ We recall that, accord- 
ing to Refs. I2H and [3(1 power-law localization induces an AT 
in any dimension provided that the decay exponent a matches 
the dimensionality d of the space. If a < (>) d, eigenstates 
tend to be delocalized (localized) in the thermodynamic limit 
but the degree of convergence is slow.— 

Based on the above facts we claim that the transition to lo- 
calization close to the origin observed in the Hamiltonian (fl} 
has its origin in the power-law localization of the eigenstates. 
We remark that this mechanism should be at work for any 
system with chiral symmetry and disorder such that a mobil- 
ity edge appears close to the origin. The Hamiltonian (fl~|i was 
chosen for the sake of simplicity. It is just one of the simplest 
representatives of the universality class associated with the AT 
in systems with chiral symmetry. 

Moreover, for a given spectral window close to the ori- 
gin, we expect the exponent a describing the decay to in- 
crease with the disorder strength W. The AT will occur for 
a W — W c so that a ~ d = 3. For stronger disorder 
W > W c eigenstates must be eventually exponentially lo- 
calized, though they may still possess a power-law tail for dis- 
tances smaller than the localization length. 

We have tested this conjecture by analyzing the typical de- 
cay of eigenstates in the window 10~ 4 < e < 10 -3 for 
L = 20 and different Ws (similar results are obtained if W 
is kept fixed and the energy window is modified) through the 
following correlation function;^ 

g(r) = <Mfe,J/hfc|)) . r - \r, - r 3max | , (4) 
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FIG. 3: The correlation function g(r) as a function of r (dashed lines) 
for a system size L — 20 and W = 12, 14, 16, 18, 20, and 22 from 
bottom to top. The numerical curves were fitted to a curve of the 
form g(r) = \n(Ar a + B) (solid lines). The best-fitting parameter 
a is represented as a function of W in the inset. In all cases the 
statistical error (not shown) is ~ a/50. 



where j max denotes the lattice site with the largest amplitude 
and the angular brackets stand for the ensemble average. In 
order to examine the asymptotic decay of g(r) (we restrict 
ourselves to the range 8 < r < 22), we fit the numerical re- 
sults to a curve of the form g(r) — ln(Ar a + B) with A, B, 
and a fitting parameters. The parameters A and B describe 
finite-size effects and remnants of the eigenstate core (nonde- 
caying) part. The exponent a (depicted in the inset of Fig. 3) 
provides us with valuable information about the decay of the 
eigenstates. 

The excellent agreement (see Fig. 3) between the numeri- 
cal results and the fitting curve indicates that eigenstates are 
power-law localized. Moreover, the exponent a controlling 



the decay increases with W and tends to the conjectured value 
a = 3 as we get close to the critical disorder W c ~ 18. 
Thus chiral symmetry combined with strong disorder induces 
power-law localization and eventually an AT close to the ori- 
gin. This is the most relevant result of the paper. 

We remark that the relation between chiral symmetry and 
power-law localization has already been established in the 
context of QCD (see Ref. 18 and references therein). Specif- 
ically, it was found that the low-lying eigenstates of the QCD 
operator (with gauge configurations given by the instanton 
liquid model) are power-law localized due to the long-range 
behavior of certain nonperturbative solutions (instantons) of 
the classical equations of motion. We speculate that a similar 
mechanism is at work in our case. 

A final comment is in order: the claim that the power-law 
localization is the precursor of the AT is not in contradiction 
with the multifractal structure of the eigenstates. Since g(r) 
is averaged over many realization of disorder, all multifrac- 
tal fluctuations are washed out. What remains is the smooth 
skeleton of the eigenstate which, in our case, has a power-law 
tail. 

To conclude, we have investigated the interplay between 
chiral symmetry and Anderson localization in a 3D Anderson 
model with off-diagonal disorder. Close to the origin we have 
found a mobility edge induced by the chiral symmetry with 
properties quantitatively different from those of a standard 3D 
AT. This chiral AT can be traced back to the power-law nature 
of eigenstates close to the origin. Our results are relevant 
to chiral systems undergoing a metal-insulator transition as 
a function of disorder. Typical examples include the QCD 
Dirac operator around the chiral phase transition and highly 
doped semiconductors where off-diagonal bond disorder 
plays a dominant role. 
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